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ABSTRACT 



The decay constants of pseudoscalar heavy mesons of 0~ state are computed by means of the 
relativistic (instantaneous) Salpeter equation. We solved the full Salpeter equation without 
making any further approximation, such as ignoring the small component wave function. 
Therefore, our results for the decay constants include the complete relativistic contributions 
from the light and the heavy quarks. We obtain F Ds « 248 ± 27 , Fjj ~ 230 ± 25 (D°, D±), 
F Bs ps 216 ± 32 , Fg ~ 196 ± 29 (S°, J B ± ), F Bc « 322 ± 42 and F Vc w 292 ± 25 MeV. 
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1 Introduction 



The decay constants of mesons are very important quantities. The study of the decay constants has 
become an interesting topic in recent years, since they provide a direct source of information on the 
Cabibbo-Kobayashi-Maskawa matrix elements. In the leptonic or nonleptonic weak decays of B or D 
mesons, the decay constants play an important role. Further, the decay constant plays an essential role 
in the neutral D — D or B — B mixing process. 

Up until now, the only experimentally obtained values of the decay constants are those of Fjj+ 
and F Ds . The first value is F D + = 300±i|g±|g MeV by BES [1], with very large uncertainties. The 
experimental values of Fd 3 have been obtained from both D s — > /j,u^ and D s — > tv t branching fractions 
by many experimental collaborations (Refs. [2, 3, 4, 5, 6, 7, 8, 9, 10, 11]). They are shown in Table 1. 
The central values from various experiments range from 194 to 430 MeV. The experimental uncertainties 
in each experiment are large, even in the most recent measurement, by ALEPH [11] (Fd s = 285 ± 19 ±40 
MeV), which has the smallest uncertainty. Further, also in ALEPH's measurement, the contribution from 
the decay D s — > fiv^j is ignored. Unlike the tree level case which is helicity-suppressed, this radiative 
decay does not have the helicity suppression. Therefore, this radiative decay may contribute several per 
cent to the branching ratio [12], and may thus cause a sizeable change in the value of the decay constant 
Fd b . Fortunately, new experiments such as Belle, BaBar, Tevatron Run II and CLEO-c will give us a 
wealth of precision data for B and D mesons soon, and will determine the decay constants to a few per 
cent. 

Many theoretical groups are working on the calculation of the decay constants, using different models, 
for example, lattice QCD, QCD sum rules, and the potential model. In Fig. 1 (taken from Ref. [13]), 
as an example, the world average of the quenched lattice results for Fb b [14, 15] is shown. From this 
Figure, we conclude that they give a stable estimate of the decay constant Fb s over several years, but the 
uncertainties are not small and still remain unchanged over the last several years of work. More precise 
predictions are still not available. Decay constants of other mesons calculated by lattice methods face 
the same problem as Fb 3 , the uncertainties are still large. Precise experimental results with uncertainties 
of only a few per cent will be obtained soon. Therefore, advancing the lattice QCD calculations is an 
urgent task in the future. Calculations with different models are and will continue to be needed, as a 
means to cross-check the lattice results. 
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Table 1: Summary of the experimental determinations of the decay constant Fd s 



Ref. 


F Da (MeV) 


[2] WA75 '93 


232±45±52 


[3] CLEO I '94 


344±37±52±42 


[4] E653 '96 


194±35±20±14 


[5] L3 '97 


309±58±33±38 


[6] DELPHI '97 


330±95 


[7] BES '98 


430tJ^±40 


[8] CLEO II '98 


280±19±28±34 


[9] BEATRICE '00 


323±44±12±34 


[10] OPAL '01 


286±44±41 


[11] ALEPH '02 


285±19±40 
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Figure 1: World average of the quenched lattice estimates of Fb s 



In this letter, we present results of a relativistic calculation of decay constants in the framework of 
full Salpeter equation. The full Salpeter equation is a relativistic equation describing a bound state. 
Since this method has a very solid basis in quantum field theory, it is very good in describing a bound 
state which is a relativistic system. In a previous paper [16], we solved the instantaneous Bethe-Salpeter 
equation [17], which is also called full Salpeter equation [18]. After we solved the full Salpeter equation, 
we obtained the relativistic wave function of the bound state. We used this wave function to calculate 
the average kinetic energy of the heavy quark inside a heavy meson in 0~ state, and obtained values 
which agree very well with recent experiments. We also found there that the relativistic corrections are 
quite large and cannot be ignored [16]. In this letter we use this method to predict the values of decay 
constants of heavy mesons in _ state. 



2 



2 Decay constants of State 

In this Section, we will calculate the decay constants of heavy mesons in 0~ state by using the full Salpeter 
method. In the previous paper [16], we wrote the relativistic wave function of CT state as: 



/Vi(g) + M H <p 2 (q) - i^Ml — ~, + irfvM 



7 5 • (1) 



{m q uJ Q + m Q uj q ) {m q uJ Q + m Q u q ) 

where m Q , m q and Mh are the masses of heavy quark, light quark, and the corresponding heavy meson, 
respectively; p Q and p q are the momenta of the constituent quarks, and P the total momentum of the 
heavy meson, q is the relative momentum of the meson defined as 

q = p Q -a 1 P = a 2 P-p q , 



where 



the ujq and uj q are defined as 



m Q 

on = ; , « 2 = 



m Q +m q ' m Q + m q 



U Q = \l m Q ~ 9±' w <? = y m q ~ ih 



where the orthogonal part q± of momentum q is defined as 

?" = gf + <£, 

gf = (P ■ q/MfoP" , <£ = if - gf . 

In the center-of-mass system of the heavy meson, gii and q± turn out to be the usual components (qo,0) 
and (0, q), and ujq = (m 2 Q + q 2 ) 1 / 2 and uj q = (m 2 + q 2 ) 1 / 2 . Wave functions f\(q) and <p2(q) represent 
the eigenfunction of the heavy meson obtained by solving the full Salpeter equation. They will fulfill the 
normalization condition 



4y>! 5)^ g)M£ -9 L + -9 i + " g o ^ = 2M H . (2) 



Several input parameters are needed when we solve the full Salpeter equation. They are similar to 
those in the usual potential model. In Ref. [16], we obtained the following best fit values of the input 
parameters by fitting the mass spectra for heavy mesons of 0~ states: 

a = e = 2.7183, a = 0.06 GeV, V = -0.60 GeV, A = 0.2 GeV 2 , A QC d = 0.26 GeV and 
m b = 5.224 GeV, m c = 1.7553 GeV, m s = 0.487 GeV, m d = 0.311 GeV, m u = 0.305 GeV. 



With this parameter set, we solved the full Salpeter equation and obtained the eigenvalues and the 
eigenfunction of the ground heavy 0~ states. We will not show here how the full Salpeter equation is 
solved and what the calculated mass spectra are, interested reader can find them in Ref. [16]. We can 
use the obtained eigenfunction of heavy mesons to calculate the decay constant Fp. The decay constant 
is defined as 

(0\qiWq2\P) = iFpPp, (3) 
which can be written in the language of the Salpeter wave functions as: 



(0\q W Q\P) = %VN c jTr [ 7m (1- 75)^5o®7^3] =^VN C P„ J t£^M0)- (4) 



dq 



(27T) 



(27T) 

Therefore, we have 

dq 



Fp = ^cI j^Mq), (5) 



and the calculated values of decay constants are displayed in Table 2. 

In Table 3, we show the theoretical uncertainties of our results for the decay constants. These 
uncertainties are obtained by varying all the input parameters simultaneously within 10% of the central 
values, and taking the largest variation of the decay constant. 

Table 2: Decay constants of heavy 0~ meson (in GeV) as predicted by the relativistic Salpeter method 





Fb s 


F Bd 


f b u 


F Vc 


Fd s 


F Dd 


f d u 


322 


216 


197 


196 


292 


248 


230 


230 



Table 3: The theoretical relative uncertainties, obtained as explained in the text, in per cents (%). 





B c 


B s 


B d 


B u 


Vc 


D s 


D d 


D u 


AFp/Fp 


±13 


±15 


±15 


±15 


±8.6 


±11 


±11 


±11 



In Table 4, for comparison, we show recent theoretical predictions for the decay constants as obtained 
by other methods. For example, we display the recent values from relativistic potential model (PM) [19] 
based on the quasi-potential approach; most recent value of Fp from another version of using Bethe- 
Salpeter method (BS) [20], which is also a relativistic result; recent values from the averaged lattice 
results both in quenched (AQL) and unquenched (AUL) approximation [15]; most recent values from 
quenched lattice (QL) QCD [21, 22] and unquenched lattice (UL) QCD [23]; and values from QCD Sum 
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Table 4: Recent calculations by other methods. Here PM means Potential Model, BS means Bethe-Salpeter 
method, QL means Quenched Lattice, AQL means Average Quenched Lattice, UL means Unquenched Lattice, 
AUL means Averaged Unquenched Lattice, QSR means QCD Sum Rules. In Ref. [21], the uncertainties are 
statistical, systematic within the Nf = 2 partially quenched approximation, the systematic errors due to partial 
quenching and the missing virtual strange quark, and an estimate of the effect of chiral logarithms, respectively. 
In Ref. [23], the uncertainties are from statistics, chiral extrapolation and systematics. 



Ref. 



F B , or F Bu 



F D , or Fn 



PM[19] 
BS[20] 



196±20 



178±15 
192 



266±25 



243±25 



QL[21] 217(6)(±|)(± 9 3 )(±i 7 ) 
QL[22] 

AQL[15] 200±20 
UL[23] 

AUL[15] 230±30 



190(7)(i-)(i-)(i«) 

173±23 
190(14) (07) (19) 
198±30 



252±9 
230±14 

250±30 



215(6)(^)(l« 3 )(+ 4 ) 



203±14 



226±15 



QSR[24] 
QSR[25] 



236±30 



203±23 
206±20 



235±24 



204±20 
195±20 



Rules (QSR) [24, 25]. As can be seen from Tables 2 and 4, our values of the decay constants by solving 
the Salpeter equation, agree with these recent results by other methods. In particular, they agree very 
well with the recent average of the unquenched lattice QCD (AUL) [15]. Our value Fq s « 248 GeV is 
smaller than the most recent experimental central value, the ALEPH's value Fd b ~ 285 ± 19 ± 40, but 
still within the experimental uncertainties. 



Table 5: Ratios Fs s /FB d , Fr>JFo d and the Grinstcin ratio by this work and by other methods. In Ref. [26], the 
first and second uncertainty arc the statistical and the systematic errors. 



Ref. 



IK 



~Fd7 



Rl 



This work 
PM[19] 



1.10±0.01 
1.10±0.21 



1.08±0.01 
1.09±0.22 



1.02±0.02 
1.01±0.40 



QL[21] 
UL[26] 
AQL [15] 
AUL [15] 



1.16(1)(2)(2)(1*; 

1.15±0.03 
1.16±0.05 



1.14(l)(il)(3)(l) 

1.12±0.02 
1.12±0.04 



1.02(2)(4)(4)(1*) 
1.018±0.006±0.010 
1.03±0.05 
1.04±0.08 



QSR[27] 



1.16±0.05 



1.15±0.04 



1.01±0.08 



There are other interesting quantities such as the ratios of decay constant Fb s / Fs d , Fn s / Fr> d , and 
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the Grinstein ratio [28] defined as 

Fn Fn 

m = (-^)/(-^), (6) 

which is a quantity sensitive to the light quark mass m s . In Table 5 we show our values of these 
ratios and some values obtained by other methods in recent literature. Our uncertainties come from the 
aforementioned 10 per cent changes of the parameters. The uncertainties of the ratios of decay constants 
of Ref. [19] are large. This is so because the authors of Ref. [19] did not give the uncertainties for these 
ratios. We estimated the uncertainties of these ratios on the basis of their given uncertainties of the decay 
constants. In the same way we estimated the uncertainties of the Grinstein ratio of other references shown 
in Table 5, with the exception of those of Ref. [26]. From Table 5 one can see that our values of ratios 
FB 3 /FB d and Fr,jFr, d agree with these recent theoretical results. In particular, our central values are 
very close to those of the relativistic potential model [19], and our central value of the Grinstein ratio 
Rl = 1.02 agrees well with the values estimated by other methods. 

In conclusion, we calculated the decay constants of heavy 0" mesons by means of the relativistic 
Salpeter method. We obtained F Da » 248 ± 27 , Fp ~ 230 ± 25 (L>°,L> ± ), F Bs » 216 ± 32 , 
F B ~ 196 ± 29 ( J B°,S ± ), F Bc ~ 322 ± 42 and F Vc ss 292 ± 25 MeV. 

The work of C.S.K. was supported in part by CHEP-SRC Program, in part by Grant No. R02-2003- 
000-10050-0 from BRP of the KOSEF. The work of GW. was supported by BK21 Program. The work 
of G.C. was supported by FONDECYT (Chile) grant no. 1010094. 

References 

[1] J. Z. Bai et al., BEPC BES Collaboration, Phys. Lett. B429 (1998) 188. 
[2] S. Aoki et al., WA75 Collaboration, Prog. Theor. Phys. 89 (1993) 131. 
[3] D. Acosta et al., CLEO Collaboration, Phys. Rev. D49 (1994) 5690. 
[4] K. Kodama et al., E653 Collaboration, Phys. Lett. B382 (1996) 299. 
[5] M. Acciarri et al., L3 Collaboration, Phys. Lett. B396 (1997) 327. 



[6] DELPHI Collab., DELPHI note 97-105 CONF 87 (1997). 

[7] J. Z. Bai et al., BEPC BES Collaboration, Phys. Rev. D57 (1998) 28. 

[8] M. Chada et al., CLEO Collaboration, Phys. Rev. D58 (1998) 032002. 

[9] Y. Alexandrov et al., BEATRICE Collaboration, Phys. Lett. B478 (2000) 31. 

[10] G. Abbiendi et al., OPAL Collaboration, Phys. Lett. B516 (2001) 236. 

[11] A. Heister et al., ALEPH Collaboration, Phys. Lett. B528 (2002) 1. 

[12] Guo-Li Wang, Chao-Hsi Chang and Tai-Fu Feng, hep-ph/0102251; hep-ph/0102293. 

[13] Damir Becirevic, hep-ph/0310072. 

[14] J. Flynn, Nucl. Phys. Proc. Suppl. 53 (1997) 168; T. Onogi, ibid 63 (1998) 59; T. Draper, ibid 73 
(1999) 43; S. Hashimoto, ibid 83 (2000) 3; C. W. Bernard, ibid 94 (2001) 159; N. Yamada, ibid 119 
(2003) 93. 

[15] S. Ryan, Nucl. Phys. Pro. Suppl. bf 106 (2002) 86. 

[16] C. S. Kim and Guo-Li Wang, Phys. Lett. B584 (2004) 285. 

[17] E. E. Salpeter and H. A. Bethe, Phys. Rev. 84, (1951) 1232. 

[18] E. E. Salpeter, Phys. Rev. 87, (1952) 328. 

[19] D. Ebert, R. N. Faustov and V. O. Galkin, Mod. Phys. Letts. A17 (2002) 803. 

[20] Zhi-Gang Wang, Wei-Min Yang and Shao-Long Wan, Phys. Lett. B584 (2004) 71. 

[21] C. Bernard et al., Phys. Rev. D66 (2002) 094501. 

[22] A. Juttner and J. Rolf, Phys. Letts. B560 (2003) 59. 

[23] N. Yamada et al., Nucl. Phys. Pro. Suppl. 106 (2002) 397. 

[24] Stephan Narison, Phys. Letts. B520 (2001) 115. 

[25] A. A. Penin and M. Steinhauser, Phys. Rev. D65 (2002) 054006. 

7 



[26] T.Onogi et al., Nucl. Phys. Proc. Suppl. 119 (2003) 610. 
[27] Stephan Narison, Phys. Letts. B322 (1994) 247. 
[28] B. Grinstein, Phys. Rev. Lett. 71 (1993) 3067. 



8 



